Diversity of flux avalanche patterns in superconducting films 
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The variety of morphologies in flux patterns created by thermomagnetic dendritic avalanches in 
type-II superconducting films is investigated using numerical simulations. The avalanches are trig- 
gered by introducing a hot spot at the edge of a strip-shaped sample, which is initially prepared 
in a partially penetrated Bean critical state by slowly ramping the transversely applied magnetic 
field. The simulation scheme is based on a model accounting for the nonlinear and nonlocal electro- 
dynamics of superconductors in the transverse geometry. By systematically varying the parameters 
representing the Joule heating, heat conduction in the film, and heat transfer to the substrate, a 
wide variety of avalanche patterns is formed, and quantitative characterization of areal extension, 
branch width etc. is made. The results show that branching is suppressed by the lateral heat diffu- 
sion, while large Joule heating gives many branches, and heat removal into the substrate limits the 
areal size. The morphology shows significant dependence also on the initial flux penetration depth. 

PACS numbers: 74.25.Ha, 68.60.Dv, 74.78.-w 



I. INTRODUCTION 

The gradual flux penetration in type-II superconduct- 
ing films is occasionally interrupted when large amounts 
of magnetic flux rush in from the edges, forming complex 
dendritic patterns inside the specimen. Experiments us- 
ing magneto-optical imaging have revealed that such den- 
dritic avalanches takeplace in numerous materials, e.g., 
NbEl YBa 2 Cu 3 7 _,PEl MgB 2 pE Nb 3 Sn,EI YNi 2 B 2 CP 
PlP and NbNEHl see also Ref. [11] for a review. The simi- 
larity in the avalanche morphology in so different mate- 
rials strongly suggests that the origin of the phenomenon 
is of a universal nature. Indeed, it is today widely agreed 
that the dendritic instability originates from a thermo- 
magnetic breakdown in the superconductor This can 
occur when a small temperature fluctuation locally weak- 
ens the pinning of the vortices, causing some magnetic 
flux to advance into the superconductor. Since moving 
flux releases heat, the local temperature then increases 
further, and a positive feedback loop is formed, which can 
lead to a rapid runaway in the temperature and magnetic 
flux propagation. 

In superconducting films placed in perpendicular mag- 
netic fields, the nonlocal electrodynamics complicates the 
theoretical description of the instability. Thus, an ana- 
lytical treatment was largely delayed compared to the 
bulk case, and in particular, the conditions for insta- 
bility onset, and the origin of a fingering nature of the 
avalanches in the film case were explained only relatively 
recently.^ES] To follow the complete time evolution of 
dendritic avalanches, including the cascades of branching 
events, only numerical simulations have proved successful 
and indeed produced very realistic result sP^Hm However, 
so far little work has been done to systematically charac- 
terize these patterns, and to our knowledge, no effort was 
made trying to identify how various physical parameters 
influence the morphology of the avalanches. 



Such characterization is not readily done experimen- 
tally since in practice it is difficult to vary the material 
parameters independently. Hence, the most feasible way 
to make a systematic investigation is to carry out a nu- 
merical simulation study. In this paper we present results 
of such a systematic study of the morphology of den- 
dritic flux patterns resulting from the thermomagnetic 
instability in superconducting films. First, we identify 
the dimensionless parameters that enter the governing 
equations. Then, avalanches are nucleated at an edge, 
where partially penetrated critical-states serve as initial 
conditions, and we follow the avalanche evolution until it 
ends in a frozen flux pattern. The final patterns created 
using different parameter values are compared, analyzed 
quantitatively and discussed. 

The paper is organized as follows: Section II describes 
the model and introduces the key dimensionless param- 
eters. Section III reports and discusses the results of our 
simulations where those parameters are systematically 
varied. In addition, the dependence of the avalanche 
patters on the initial critical-state is investigated. The 
results are summarized in Sec. IV. 



II. MODEL 

Consider a long superconducting strip of half- width w 
and thickness d <C w, in thermal contact with a sub- 
strate, see Fig. [T] In the numerical simulations the fol- 
lowing strategy is used: Starting from a zero-field-cooled 
strip, we slowly ramp (increase) the applied magnetic 
field H a , while the thermal feedback is turned off. This 
ensures that the flux penetration is gradual and the tem- 
perature is everywhere equal to the substrate temper- 
ature To. The spatial distributions of flux density B z 
and sheet current J develop in agreement with the crit- 
ical state model for the transverse geometry.^ The field 
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FIG. 1. Sketch of the sample configuration. A long super- 
conducting strip of thickness d and half- width w, is in thermal 
contact with a substrate. Over a region of width I near both 
edges a critical state with sheet current J — J c and flux den- 
sity B z exists as a hot-spot (red dot) is introduced to trigger 
an avalanche. 



ramp is stopped when the flux penetration reaches the 
depth, Z, where l/w ~ 0.1, and that state is taken as ini- 
tial condition for the upcoming avalanche. The thermal 
feedback is then turned on, and a heat pulse is applied 
to a small region near the edge, thus triggering the in- 
stability. The same protocol is repeated using different 
sample parameters, and for each run we analyze the size 
and morphology of the flux distribution frozen in the su- 
perconductor after the avalanche. 

Before describing the physical model in detail, let us 
introduce the units that will be used. Time is expressed 
in units of to = fiodw/ p n , where /io is the vacuum mag- 
netic permeability and p n is the normal-state resistivity 
of the superconductor. The sheet current J and mag- 
netic field B z //jLq are both expressed in units of the zero- 
temperature critical sheet current, J c o = dj c o. The tem- 
perature T is normalized to the critical value T c , and the 
spatial coordinates (x, y) are in units of w. Consequently, 
the Faraday law, B = —V x E, implies that the electric 
field E is normalized by p n jc0- 

We consider a superconductor described by a non- 
retarded relationship between the electrical field and the 
current density as follows 



E = pJ, p 




if J > J c or T > 1, 
otherwise . 



(i) 



Here J = | J| while n is the creep exponent. The temper- 
ature dependencies of n and J c are taken aP^ 



J C (T) = 1-T, 
n(T) = rti/T + n . 



(2) 



The thermal properties of the superconducting film are 
specified by the temperature dependence of the specific 
heat c, thermal conductivity ft, and the coefficient of 
heat removal to the substrate h. We use here an ap- 
proximation applied successfully in previous modelling 



work, namely that all three parameters have cubic T- 
dependenciesPSl 

Based on this, the evolution of the local temperature 
is described by the dimensionless diffusion equation, 



t = aS7 2 T - (3(T - T ) + <yT~ 3 J • E , 



(3) 



where To is the normalized substrate temperature. The 
coefficients 



dp K p h ,MoJ? 
a= , fi = w , 7 = wd— — , 

W p n C p n C C T c 



(4) 



are also dimensionless, and involve c, ft and h evaluated 
at T c . The parameter a has the meaning of a normal- 
ized coefficient of thermal diffusion, and characterizes the 
smearing of T due to heat conduction within the film. 
The f3 characterizes the heat flow to the substrate, while 
7 measures the positive feedback due to Joule heating. 
For 7 = the evolution of the magnetic flux and temper- 
ature distributions are decoupled and the film is always 
thermomagnetically stable. 

Note that in this representation the parameter space 
is significantly reduced compared to the dimensional de- 
scription, since all the ft, /i, c, T c , j c0 , and p n have been 
combined into three dimensionless parameters. 

To calculate the electromagnetic behavior we express 
the sheet current through the local magnetization, g, as 



J x = dg/dy, J y 



-dg/dx. 



(5) 



Outside the sample, g vanishes by definition. In an infi- 
nite or periodic space, the Biot-Savart law for a thin film 
has a simple expression in the Fourier space. Inverting it 
and taking the time derivative yields 16 



9 = ?- 




(6) 



where T and T 1 are the forward and inverse Fourier 



transform, respectively, and k = y k% + k% is the in-plane 

wave- vector. Inside the sample, the right-hand side of 
Eq. Q is found from the Faraday law in combination 
with the material law, Eq. 0. Outside, B z is found 
implicitly by an iterative scheme requiring that g = 0. 
More details about the simulation procedure are found 
in Ref . [16j 

The sample occupies a 2 x 2 square, with periodic 
boundary conditions in one direction. The other direc- 
tion is extended with extra space in order to take into 
account the magnetic field outside the strip. The total 
area of the simulation covers a rectangle of size 3x2, 
which is discretized on a 768 x 512 equidistant grid. 

The simulations were carried out using To = 0.15, and 
the creep exponent was set to n = 20/T — 10. The ini- 
tial state was prepared with a ramp rate of H a — 10 -5 , 
and the hot spot triggering the avalanche had a size of 
0.01 x 0.01 (a cross consisting of 5 grid points) and a 
temperature of T = 1.5. During the avalanche the field 
ramp was stopped. 
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FIG. 2. Variety of avalanche flux distributions B z arranged in a matrix of panels with constant a in the rows, and with 
constant 7 in the columns. The image brightness represents the magnitude of B z , as in experimental magneto-optical images, 
and the strip edges are seen as a pair of bright horizontal lines while the dark area between them is the flux-free Meissner state 
region of the superconducting strip. 
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III. RESULTS AND DISCUSSION 

The parameter space in these simulations is three- 
dimensional, and we present and discuss first the various 
types of flux patterns produced with j3 = 0.1, i.e., keep- 
ing a fixed coefficient of heat transfer to the substrate. 
The initial flux penetration depth was / = 0.11 when a 
hot spot was introduced to trigger an avalanche. 



A. Varying a and 7 

Figure 2 displays a collection of final-state flux 
avalanche patterns obtained using values for a and 7 in 
the range 10 -5 — 10 -2 and 1 — 20, respectively. Evidently, 
within this part of parameter space one finds patterns 
which show striking similarities with those observed in 
experiments In addition, we find some type of pat- 
terns not previously reported. 

The figure shows that with an increasing Joule heating, 
7, and a constant a, the avalanches quickly become larger 
in size and get more and more branches. At the same 
time, the branch width remains essentially the same, al- 
though for increasing 7 the main trunk near the edge gets 
steadily wider. 

When instead keeping 7 constant and increasing the 
in-plane heat diffusion, a, one sees that the number of 
branches decreases significantly. At the same time, the 
overall avalanche size is not varying much. At maximum 
value, a = 10 -2 , the avalanches can hardly be character- 
ized as dendritic, but rather as a soft protrusion. Such 
non-dendritic shapes is a consequence of the thermal dif- 
fusion being almost as fast as the electromagnetic prop- 
agation. 

The simulated I^-maps reveal also several other in- 
teresting features. One is that one never finds branch- 
ing to take place inside the critical-state region near the 
edge. Another is that even in the largest avalanches, 
the branches rarely propagate past the sample center. 
Furthermore, branches never overlap, thus appearing to 
repel each other. Yet another is that the inner part of 
the branches most often has the highest flux density, but 
in some cases they show a dark low-I^ core. Note also 
that the flux density is always reduced in the avalanche 
trunk as well as outside the sample, next to the root lo- 
cation. All these features correspond very well with ex- 
perimental observations using magneto-optical imaging, 
thus demonstrating a detailed correspondence between 
our simulations and reality. 

A feature seen in Fig. [2] which so far was not reported 
neither from experiments nor from previous simulations, 
is the appearance of quasi-periodic side branches in the 
panels with small a. The origin of these branches is not 
yet clear, but a similar phenomenon is known from other 
systems, e.g., formation of dendrites in crystal growth.^ 

The panels in the upper left corner of Fig. |2j corre- 
sponding to small 7 and large a, do not show any visible 
trace of avalanche activity. This indicates the existence 



of a region in parameter space where the system is stable 
towards even large perturbations, such as local heating 
above T c . Hence, a guiding line for design of, e.g., su- 
perconducting power-devices is to make materials and 
dimensions such that 7 is small and a is large. Note that 
the transition between stable and unstable behavior is a 
crossover rather than a sharp phase boundary. 

Finally, note also that the flux structures in Fig. [2] 
are more symmetric and with straighter branches than 
the avalanche structures reported fr om p revious simula- 
tion work using the same formalism! 16 * 17 * This difference 
is due to the fact that in previous work static disorder 
was introduced in the critical current density whereas 
here the sample was perfectly uniform. Regarding the 
role of disorder we therefore conclude that (i) Branching 
flux patterns may appear even in uniform materials, in 
agreement with theoretical predictions of the fingering 
instability in a uniform systemJ^^ (ii) It is the non- 
linear and nonlocal electrodynamics, rather than spatial 
disorder, that is responsible for the branching of the flux 
structures, (iii) The presence of disorder affects the exact 
path along which the hot branches propagate. Straight 
main branches with quasi-periodic side branches in a uni- 
form sample become wiggly and lose periodicity in a dis- 
ordered sample, (iv) Figure [2] shows no sign of symmetric 
bifurcations, suggesting that the process of symmetric tip 
splitting found in Ref. 17 is also driven by disorder. 

As a quantitative measure of the effect of varying 7 
and a, Fig. [3] shows the total avalanche size, taken as 
the area of the Meissner state region invaded by flux, 
plotted as a function of 7 for various a. In spite that 
a spans 4 decades the graphs are not very different and 
increase almost linearly with the same slope. Evidently, 
the Joule heating parameter 7 is a crucially important 




FIG. 3. Avalanche size in units of w 2 as a function of 7 
for different a. The data are extracted from the flux patterns 
displayed in the panels of Fig. [2] 
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FIG. 4. Average branch width in units of w as function of 
a for various 7, extracted from the flux patterns shown in 
Fig.0 



that the width of the individual branches in the den- 
dritic structure to a large extent is controlled by thermal 
diffusion. 



B. Varying f3 

Consider next how the avalanche morphology is influ- 
enced by varying the heat transfer to the substrate char- 
acterized by the parameter /3. Figure [5] illustrates the 
variety of flux patterns obtained for (3 ranging from 0.1 
to 100, while keeping constant the parameters a = 10 -4 , 
7 = 10, and I = 0.11. It is evident that an increase in 
P leads to avalanches of smaller overall size and thinner 
branches, and to some extent, also less degree of branch- 
ing. The behavior can be understood from the role of j3 in 
the heat diffusion equation, Eq. ([3]). During an avalanche 
the term proportional to (3 represents a heat-sink that 
stabilizes the runaway and thus tend to terminate the 
avalanche P^l A large heat sink effect is also expected to 
give thinner branches. 



for the impact of an avalanche. 

Presented in Fig. [4] is the average branch width plotted 
as a function of a for various 7. The width is evaluated 
by dividing the area of the avalanche by the length of 
the contour of B z = around the structure. The fig- 
ure shows that the width grows rapidly with increasing 
a, and is quite weakly depending on 7. This suggests 




FIG. 5. Avalanche flux patterns B z obtained for various 
/3, keeping constant a = 10 -4 , 7 = 10 and / = 0.11. These 
panels show only the lower half of the strip. 



C. Varying / 

Finally, we consider how the avalanches depend on the 
penetration depth, Z, of the critical-state background flux 
distribution. Shown in Fig. [6] are four images of the strip 
where I varies from 0.04 to 0.3 while a, f3 and 7 are kept 
constant. One sees that large I also gives large avalanches 
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FIG. 6. Flux distributions B z for avalanches nucleated at 
different critical-state penetration depths / = 0.04, 0.1, 0.2, 
and 0.3, with a 10" 4 , f3 0.1 and 7 10. 
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and many branches. This results is in agreement with 
previous experiment by Bolz et alP 

The images seen in Fig. [6] show also a striking resem- 
blance with previous experimental and sim ulation results 
at different substrate temperatures To The results of 
this work suggests that the increasing avalanche size with 
Tq is not because the avalanche propagation is sensitive 
to To, but rather because it is sensitive to the flux pen- 
etration depth prior to the avalanches, which increases 
with T . 

IV. SUMMARY 

In this work we have used numerical simulations to in- 
vestigate how the dimensionless material parameters af- 
fect the size and morphologies of dendritic flux avalanches 
in superconducting films. The avalanches were nucleated 
at the edge of a strip by a heat pulse and by varying 
the parameters we have reproduced a wide range of mor- 
phologies previously seen experimentally. We find that 
increasing the normalized coefficient of heat diffusion a 
gives fewer and wider branches while strong heat removal 



to the substrate, quantified by /?, cause smaller total 
avalanche size and thinner branches. Increasing values 
for the Joule heating parameter 7 gives more branches 
and larger avalanches. Parameter combinations with 
large a and small 7 proved to be stable towards large 
perturbations. Quasi-periodic side branches were seen 
in the dendritic structures with small a, and their exis- 
tence was attributed to the absence of spatial disorder. 
Finally, we found that the avalanche morphology is sen- 
sitive to the initial critical state before the avalanche, 
namely, deeper initial penetration gives larger and more 
branching avalanches. 

For future work, it will be of interest to use the present 
dimensionless description to study how the linear stabil- 
ity diagram of superconductors is affected by varying the 
parameters a, /?, and 7. 
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